A Study on the Noise Threshold of Fault-tolerant Quantum Error Correction 
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Quantum circuits implementing fault-tolerant quantum error correction (QEC) for the three qubit 
bit-flip code and five-qubit code are studied. To describe the effect of noise, we apply a model based 
on a generalized effective Hamiltonian where the system-environment interactions are taken into 
account by including stochastic fluctuating terms in the system Ifamiltonian. This noise model 
enables us to investigate the effect of noise in quantum circuits under realistic device conditions and 
avoid strong assumptions such as maximal parallelism and weak storage errors. Noise thresholds of 
the QEC codes are calculated. In addition, the effects of imprecision in projective measurements, 
collective bath, fault-tolerant repetition protocols, and level of parallelism in circuit constructions 
on the threshold values are also studied with emphasis on determining the optimal design for the 
fault-tolerant QEC circuit. These results provide insights into the fault-tolerant QEC process as 
well as useful information for designing the optimal fault-tolerant QEC circuit for particular physical 
implementation of quantum computer. 



I. INTRODUCTION 

Recent developments in the theory of quantum com- 
putation have generated significant interest in utihzing 
quantum mechanics to achieve new computational capa- 
bility |l| . A quantum computer can outperform its classi- 
cal counterpart and provide efficient ways to solve many 
important problems. However, the intrinsic sensitivity of 
a quantum superposition state to imperfect operations 
and interactions with its surrounding environment pro- 
hibits the realization of a scalable quantum computer. To 
combat the inevitable errors and decoherence of quantum 
states during the process of computation, quantum error 
correction (QEC) and fault-tolerant methods of quantum 
computation have to be applied in the construction of 
large-scale quantum computers. It has become clear that 
the future of robustly storing and manipulating quantum 
information rely upon the success of fault-tolerant quan- 
tum error correction H, H, ^ . 

Fault-tolerant methods combined with concatenated 
coding yield the threshold result, that states if the 
noise level per elementary operation is below a threshold 
value, then arbitrarily long quantum cornputation can 
be achieved using faulty components 0, IE S 111- Us- 
ing a t-error correcting code, fault-tolerant circuits con- 
structed from faulty gates with error rate e can achieve 
a logical error rate of ©(e*"*"^) per logical gate. This 
fact together with the concept of concatenated coding 
provide a method for possible large-scale quantum com- 
putation, and can lead to the realization of a scalable 
quantum computer. Therefore, it is important to study 
fault-tolerant methods and estimate the noise threshold 
values. In addition, the noise threshold indicates the tol- 
erable noise level in a certain quantum circuit, and pro- 
vides a benchmark for the efficiency of QEC circuits. 



'Electronic address: 'silbey @mit . edu] 



A number of theoretical estimates of noise threshold 
and improvements for the efficiency of QEC circuits have 
been proposed |5|, i^, i7|, iJ, 0, HO, ! however, these stud- 
ies adopt ad-hoc classical stochastic noise models that ne- 
glect device details, and make strong assumptions such as 
maximal parallelism and low noise level in storing qubits. 
Realistically, these assumptions are not usually applica- 
ble, and the power of fault-tolerant QEC under realistic 
physical conditions is still unclear. Noise threshold values 
are of little use if limitations of the physical implementa- 
tion and realistic noise sources are not considered in the 
estimation. Therefore, it is of importance to study fault- 
tolerant QEC circuits using a noise model that reflects 
realistic device conditions. 

In Ref. |T2|, we applied a phenomenological noise 
model to study the effect of noise in quantum teleporta- 
tion and controUed-NOT gate operation. Starting from a 
effective system Hamiltonian that incorporates stochas- 
tic fluctuating terms to describe the effect of system- 
environment interactions, this model can describe the 
dissipative dynamics of a many-qubit system under re- 
alistic device conditions. In this paper, the same model 
is applied to investigate the performance of fault-tolerant 
QEC circuits implementing three qubit bit-flip code and 
flve-qubit code. Relatively small codes are studied be- 
cause we perform a systematic investigation on several 
variables that can affect the performance of fault-tolerant 
QEC circuits. In section ^ we flrst present the model 
Hamiltonian we used to implement quantum gates, and 
briefly review the noise model we proposed. We then in- 
troduce the fault-tolerant QEC circuits studied in this 
work in section IIIII and show our estimates of noise 
threshold in section Hvl Finally, a systematic study on 
how factors like imperfect measurement, collective bath, 
repetition protocol, and level of parallelism affect the per- 
formance of fault-tolerant QEC is presented in section 
IVl This theoretical study will be useful for the design 
and implementation of fault-tolerant QEC circuits. We 
briefly conclude our results in section IVTI 
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II. INTERACTIONS AND NOISE MODEL 

We study the performance of fault-tolerant QEC cir- 
cuits using a microscopic model described in Ref . . In 
this model, a qubit system is described by a Hamiltonian 
with a controlled part and a time dependent stochastic 
part. The general Hamiltonian of the qubit system can 
be written as {h = 1) 



- Mn 



= Measurement 



n 



=exp{-i^(X+Z)} 



H(t) = Ho(t)+h(t), 



(1) 



where the controlled Hamiltonian Ho(t) describes the in- 
teractions between qubits, and the stochastic part h{t) 
describes the fluctuations of the interactions due to the 
coupling to the environment. During the process of quan- 
tum computation, Ho(i) is controlled to implement gate 
operations, whereas h.(t) is stochastic and results in the 
decoherence of the quantum system. 

We choose to simulate fault-tolerant QEC circuits us- 
ing a model control Hamiltonian with single-qubit X, Z 
and two-qubit ZZ interactions: 

Ho (t) = J2 ^'^ ^' +J2 (^)^^ + J2 9^J (t) Z, , (2) 

where Zi and Xi arc Pauli operators acting on the i-th 
qubit, andei(t), Ji(i), and5y (i) are controllable parame- 
ters that can be turn on and off to implement desired gate 
operations. For simplicity, all gate operations are simu- 
lated using step function pulses with field strengths set 
to 1 (uniform field strengths), and the "on-time" of each 
pulses are controlled to obtain the desired unitary trans- 
formations. Note that by doing so we adopt a dimension- 
less system in which a unit time scale Ai is defined by the 
field strength e, i.e. Ai = We consider fault-tolerant 
QEC circuits composed only of single-qubit Hadamard 
gates, two-qubit controUed-NOT and controUed-Z gates, 
plus measurement of a single qubit in the computational 
basis. All these operations can be easily implemented us- 
ing the model Hamiltonian in Eq. l^J. Figure^shows the 
gate symbols and corresponding unitary transformations 
used in our simulations. Note that we adopt the ZZ type 
two-qubit coupling in our model Hamiltonian for an illus- 
trative purpose. The real form of the inter-qubit interac- 
tion depends on the controllable interactions available for 
each individual physical implementations. Nevertheless, 
our model can handle the other types of interactions as 
well, and we expect that the model Hamiltonian we use 
here can reproduce the same general physical behavior 
as other two-qubit Hamiltonians. 

To describe the effect of noise, we consider fluctuations 
on the system Hamiltonian due to system-environment 
interactions: 

i i 

where fei(<) and SJi{t) describe the time-dependent di- 
agonal and off-diagonal fluctuations on the i-th qubit. 



-^V = exp{i^Z} 



z 
z 



= exp{-i^ZiZ2} 



FIG. 1: Quantum gate symbols used to denote unitary trans- 
formations implemented with single-qubit X, Z and two-qubit 
ZZ interactions. More complicated transformations such as 
controUed-NOT and controUed-Z gates can be trivially con- 
structed using these elementary gates. 



respectively. This corresponds to stochastic single-qubit 
phase (Z) and bit-flip {X) errors on each individual 
qubit. In addition, we consider the fluctuations as ran- 
dom Gaussian Markov processes with zero mean and S- 
function correlation times described by the following set 
of equations: 



{Se,{t)SJ,{t')) 



(5J,it)) = 0, 
70 • 6ijS{t - t'), 

0, 



(3) 



where bracket () means averaging over the stochastic 
variables, and 70 and 71 describe the strength of the diag- 
onal energy fluctuations and off-diagonal matrix element 
fluctuations, respectively. For a single-qubit system, 70 
and 71 are well-defined physical quantities, i.e. 70 and 71 
are population relaxation rate and pure dephasing rate, 
respectively 0. Note that noise strengths 70 and 71 
should be interpreted as the error rate per unit time scale 
At = 1/e, where e is the strength of the control fields. 
Also notice that we treat the correlation between differ- 
ent qubits independently, which means each qubit in the 
system is coupled to a distinct environment (bath). Later 
we will remove this constraint and examine the effect of 
a collective bath on the noise threshold value. We also 
assume that the diagonal and off-diagonal fluctuations 
are not correlated. 

For simplicity, we assume that the noise strengths are 
uniform, i.e. 70 and 71 are constants. The noise strength 
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is set to be the same on all qubits at all times, therefore, 
we do not distinguish storage and gate errors. By as- 
suming that the storage and gate errors are at the same 
level, the uniform noise assumption overestimates the er- 
rors in the system. At the same time it also avoids the 
weak storage noise assumption usually made in previous 
estimates of noise thresholds. Realistically, to perform a 
quantum gate between two distant qubits in a large-scale 
quantum circuit, multiple quantum swap gate s must be 
employed to shuffle quantum states around • Our uni- 
form noise assumption reflects the physical condition in 
this scenario. Note that more complex setups, in which 
control field and noise strengths are different for each 
individual qubits can be studied with exactly the same 
method. 

Using the method described in Ref. ^3 and Eq. QJ- 
, we can derive the equation of motion for the density 
matrix of the qubit system, and numerically propagate 
the density matrix of a system with up to twelve qubits 
(bound by the size of physical memory on a personal 
computer). This method provides an efficiently way to 
simulate quantum circuits and obtain full dynamics of 
the qubit system. 



III. FAULT-TOLERANT QEC CIRCUIT 

In this paper, we study fault-tolerant QEC circuits im- 
plementing the three qubit bit-flip code and five qubit 
code. We choose to investigate these two codes, because 
they are relatively small and allow us to perform system- 
atic studies. Previous studies on the fault-tolerant QEC 
have been mainly focused on CSS codes, especially the 
CSS [[7,1,3]] code [1111 El. Because fault-tolerant en- 
coded operations on CSS codes are easy to implement, 
CSS codes are expected to be more useful for quantum 
computation than the three qubit bit-flip code and five 
qubit code. Nevertheless, since we focus on variables 
affecting the performance of fault-tolerant QEC circuits, 
we expect that results gained in our study can be applied 
to more general codes. In this section, we introduce these 
two codes and the methods we apply to perform fault- 
tolerant QEC. 



A. Fault-tolerant QEC scheme 

We adopt the fault-tolerant QEC scheme proposed by 
DiVincenzo and Shor . This protocol utilizes cat states 
and transversal controUcd-X/Z gates to detect error syn- 
dromes. The procedure can be divided into three differ- 
ent stage: (1) ancilla preparation and verification, (2) 
syndrome detection, and (3) recovery. 

To detect syndrome fault-tolerantly, ancilla qubits 
have to be prepared in maximally entangled cat states, 
and go through a verification step to ensure that magni- 
tudes of correlated multiple-qubit errors are small. For 
example, the four-qubit cat state -^(jOOOO) -I- jllH)) is 



o> 

0> 



FIG. 2: The fault-tolerant circuit for the preparation and 
verification of the four-qubit cat state. Note that the final 
result is conditioned by the outcome of the measurement at 
the end of the circuit. If the measurement outcome is zero, 
we accept the state; otherwise, the state is discarded and the 
circuit is started over again. 



necessary for the fault-tolerant QEC of the five-qubit 
code. Figure 12 shows the circuit we used to prepare and 
verify four-qubit cat states In this circuit, an extra 
qubit is used to detect correlated X errors in the cat 
state; after the measurement, only states with measure- 
ment result equals to zero are accepted. This verifica- 
tion step ensures that a single-qubit error in the circuit 
causes at most a single-qubit error in the final cat state, 
therefore the circuit fulfills the fault-tolerant condition. 
Compared to other fault-tolerant cat state preparation 
circuits P, 0], an important feature in the circuit in 
Fig. 12 is that only a projective measurement is required 
to verify the cat state fault-tolerantly. This is possible 
because the circuit takes into account the error propaga- 
tion pattern in the preparation step. 

The ancilla cat state generated by the circuit pre- 
sented in Fig. [21is used to perform transversal controUed- 
stabilizer operation to transfer information about the er- 
rors from the data qubits to the ancilla qubits. After de- 
coding the ancilla state, projective measurement is then 
applied to obtain error syndromes. Because there are 
more gates in the circuits than the number of measure- 
ments, it is reasonable to assume that measurement has 
smaller effect on the threshold result. Therefore, we as- 
sume perfect measurement. Later we will study the effect 
of measurement errors. In addition, to ensure that we do 
not accept a wrong syndrome, we must repeat syndrome 
detection and take a majority vote. Following Shor's pro- 
tocol, the following repetition scheme is used: 

Repetition Protocol A (three majority vote): 

1. Perform the syndrome detection twice. If the same 
measurement results are obtained, the syndrome is 
accepted and data qubits are corrected. 

2. Otherwise, perform one more syndrome detection. 
If any two of the three measurement results are the 
same, the syndrome is accepted and data qubits are 
corrected. 



4 



3. If all three measurement results are different, no 
further action is taken. 

This protocol is basically a simple majority vote in three 
trials. Note that the choice of the repetition protocol is 
not unique. In fact, later we will compare protocol A 
to another protocol, and show that we can improve this 
protocol to increase the efficiency of the fault-tolerant 
QEC procedure. 

Combining all these elements, we can ensure that the 
probability of generating a two-qubit error is of order e^, 
and avoid the catastrophic propagation of errors. Dur- 
ing a quantum computation using a single-error correct- 
ing code, we can perform the fault-tolerant QEC after 
each gate operations. In consequence, single-qubit errors 
generated in earlier computation and QEC steps will be 
corrected in later QEC steps. Therefore, only two-qubit 
errors will be accumulated in a rate of order e^. As a re- 
sult, we can achieve longer computation when e is small. 



B. Three qubit bit-flip code 

The three qubit bit-flip code encodes a logical qubit in 
three physical qubits using the following logical states: 



|0l> 
|1l> 



= 1000), 
= nil). 



This code is a stabilizer code with two stabilizer opera- 
tors gi = ZZI and (72 = IZZ. The three qubit bit-flip 
code corrects single bit-flip error on any of the three data 
qubits. This code does not correct phase errors, there- 
fore it is only useful when the degradation of the quantum 
state is dominated by bit-flip errors. However, we believe 
insights gained by studying this code can be applied to 
more general quantum error-correcting codes. 

Figure O shows the syndrome detection circuit for the 
three qubit bit-flip code. In this circuit, physical qubits 
are depicted by horizontal solid lines, and quantum gates 
are represented by boxes. The quantum circuit includes 
three data qubits that take a encoded state as the in- 
put, and two pairs of ancilla qubits prepared in the Bell 
state -^(lOO) -I- 1 11)), which are used to measure the syn- 
dromes. We want to point out that only limited ability to 
perform operations in parallel is assumed in constructing 
this circuit. In addition, the circuit is arranged to min- 
imize error propagation from the ancilla qubits to the 
data qubits. At the end of the circuit, two measure- 
ments, Ml and M2, are performed to obtain the error 
syndrome. After the syndrome is confirmed according to 
the repetition protocol A, we then apply the correspond- 
ing recovery action to correct the detected error. Table 
m lists the syndrome and the corresponding recovery ac- 
tions for the three qubit bit-flip code. 

Because the three qubit bit-flip code only corrects bit- 
flip errors, we only consider off-diagonal fluctuations on 
each qubit when dealing with this code (70 = 0). Note 
that the circuit does not protect against Z errors, nor 



TABLE I: Measurement results and the corresponding actions 
required to correct the error in the data qubit for the three 
qubit bit-flip code. 
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FIG. 3: A circuit implementing the fault-tolerant syndrome 
detection for the three qubit bit-flip code. 



can it prevent the generation of Z errors. To access its 
performance on controlling X errors on the data qubits, 
we only study the fault-tolerant QEC procedure when 
initially the data qubits are in the logical |0l) state. This 
selection of initial state is unrealistic, but it allows us 
to avoid uncorrectable Z errors that will ruin the QEC 
procedure. 



C. The flve-qubit code 

The flve-qubit code is the smallest quantum code 
that corrects all single-qubit errors 0, 0|. A scheme 
for fault-tolerant quantum computation using flve-qubit 
code is presented by Gottesman in Ref. 3- Here we 
adopt the representation and fault-tolerant QEC circuit 
presented by DiVincenzo and Shor in Ref. 0- Their 
implementation uses a nine-qubit system with flve data 
qubits and four ancilla qubits, which utilizes four-qubit 
cat state --i=(|0000) -I- |1111)) for syndrome detection. 
Moreover, four syndromes are detected sequentially. It is 
straightforward to simulate the syndrome detection cir- 
cuit presented in their paper using our choice of model 
Hamiltonian (Eq. [TO . 

Ideally, multiple input states have to be studied to ob- 
tain averaged performance of the QEC procedure. To 
avoid such tedious computations, we use a logical qubit 
initially in the following pure state density matrix (in the 
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{|0l),|1l)} basis): 

This state provides an averaged measure for all possible 
logical states, thus should give us a reasonable estimate 
of the averaged circuit performance. 

Note that our setup simulates a minimal circuit for 
the fault-tolerant QEC using five-qubit code with limited 
physical resources. We expect such nine-qubit system can 
be realized on a liquid-state NMR quantum computer 
using available technologies. An experimental study on 
such minimal fault-tolerant QEC circuit will be an ex- 
cellent test for our noise model, and can also provide us 
invaluable information that is essential for the design of 
large-scale quantum computers. 



IV. ESTIMATE OF NOISE THRESHOLD 

To estimate the noise threshold for a logical operation, 
we simulate a computation where fault-tolerant QEC is 
performed after each logical operation on the encoded 
qubits, and compare the magnitude of logical errors to 
the magnitude of errors generated by the same operation 
on a bare physical qubit without QEC. We use the crash 
probability Pc to describe the amount of logical errors in 
an encoded state. The crash probability is defined as the 
probability of having an uncorrectable error in the data 
qubits, and can be obtained from the fidelity of the state 
after a perfect QEC process. 

We define a computational step as a logical gate fol- 
lowed by a fault-tolerant QEC step. If the same compu- 
tational step is applied on the data qubits repeatedly n 
times, we can describe the crash probability as a func- 
tion of n, i.e. Pc = Pcin)- In general, Pc{n) satisfies an 
exponential form: 



Pc{n)^-{l-e 



-2r„ 



(4) 



We can perform simulation and compute crash probabil- 
ity at each step, Pdn). By fitting our simulation result 
to the functional form in Eq. Q), we obtain the crash 

rate constant per computational step r„ — 

n— 

In addition, we also define the crash rate constant per 



unit time Tt 



dP^jt) 



dt 



= r„/T, where r is the time 

period required to complete a computational step. Note 
that the unit time scale Ai is defined by the strength of 
control fields e. At — 1/e. 

We compute noise threshold for a quantum memory, 
where repeated fault-tolerant QEC is applied on the data 
qubits to stabilize quantum information; and logical X 
gate, where a logical X gate followed by a fault-tolerant 
QEC step are applied on the data qubits. Figure^shows 
the crash rate constants as a function of noise strength 



for the three qubit bit-flip code, as well as the results 
for the five-qubit code. In Fig. ^ we clearly see that in 
the weak noise regime, the crash rate constant is propor- 
tional to the square of the noise strength, which reflects 
the power of the fault-tolerant QEC procedures. The 
noise threshold is obtained from the critical value where 
the crash rate constant for encoded computation cross 
over with the error rate of a bare physical qubit. At 
noise strength below the threshold value, the errors in 
the encoded state accumulated slower than for the bare 
physical qubit. At noise strength above the threshold 
value, the fault-tolerant QEC provides no benefit. For 
the three qubit bit-flip code, the noise threshold is about 
2 X 10^^ for quantum memory, and about 1 x 10^'' for 
the logical X gate. 

We also perform calculations on five-qubit code. The 
five-qubit code corrects all single-qubit error, so we can 
compute the threshold for different types of noise. Table 
nil summarizes threshold values for three qubit bit-flip 
code and five-qubit code. The difference in the noise 
threshold between quantum memory and logical X gate 
is mainly due to the different basis of comparison. For 
the quantum memory, we need to compare crash rate con- 
stant per unit time Fj to the decay rate of a free physical 
qubit; however, for the X gate, we need to use the crash 
rate constant per computational step r„ . The extra log- 
ical X operation has little effect on the crash rate per 
computational step because the fault-tolerant QEC cir- 
cuit is much bigger than the circuit for the logical X gate. 
This observation suggests that other encoded single-qubit 
operations and transversal encoded two-qubit operations 
should have similar threshold values. 

We summarize the assumptions we made for these cal- 
culations: (1) stochastic X and Z noises, (2) each qubit 
coupled to a distinct bath, (3) uniform noise strength, 
(4) perfect physical |0) states as initial states, (5) perfect 
instantaneous projective measurement. Clearly, the uni- 
form noise assumption that treats gate errors and storage 
errors on the same footing is responsible for the relatively 
low threshold we obtain for the five-qubit code. In the 
next section we will examine some of these assumptions 
and discuss different factors that might affect the perfor- 
mance of fault-tolerant QEC. 



EFFICIENCY OF FAULT-TOLERANT QEC 
CIRCUITS 



In this section, we study several variables that can af- 
fect the efficiency of the fault-tolerant QEC scheme. We 
perform a systematic investigation on the performance of 
quantum memories stabilized by fault-tolerant QEC and 
aim to generate a generic picture on how these variables 
change the efficiency of fault-tolerant QEC circuits. 
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TABLE II: Summary of noise threshold values. The noise strengths (70 and 71) should be interpreted as the error rate per 
unit time scale At — 1/e, where e is the strength of the control fields. 



three qubit bit-flip code five-qubit code 



quantum memory X gate quantum memory X gate 

X errors (71) 2.1 x 10~^ 1.2 x 10"^ 4.2 x 10"" 3.5 x 10"^ 

Z errors (70) - - 5.1 x 10"'' 4.3 x 10"^ 

Both X and Z errors - - 4.7 x 10"'' 3.9 x 10"^ 





10" 



10"° 10 10 10"' 

Noise Strength (y, ) 



FIG. 4: Crash rate constants as a function of the noise 
strength. We show crash rate constants for a quantum mem- 
ory using the three qubit bit-flip code (upper-left) and flve- 
qubit code (upper-right), and for a logical X gate on the three 
qubit bit-flip code (bottom-left) and on the five-qubit code 
(bottom-right). For the five-qubit code circuits, curves for 
different types of noise are presented. To show the threshold 
result, we also present curves for the error rate of a single 
physical qubit (dotted line). The noise threshold values are 
summarized in Table ITTl 

A. Effect of imperfect measurement 

We first test the effect of imperfect measurements on 
the performance of a quantum memory stabilized using 
fault-tolerant QEC. We use the following POVM (pos- 
itive operator-valued measure) to describe an imperfect 
projective measurement on a single qubit: 

Mo = (l-r;)|0)(0|+r;|l)(l|, 
Ml = (l-r;)|l)(l|+r;|0)(0|, 

where Mq (Mi) describes events in which basis state |0) 
(|1)) is measured, and rj is the probability of measurement 
error, i.e. a projection onto the wrong basis state. Figure 
|31 shows curves for the crash rate constant per unit time 



FIG. 5: Ft as a function of noise strength for the fault- 
tolerant QEC circuit using three qubit bit-flip code at differ- 
ent level of measurement errors. The error rate of a single 
physical qubit is also shown (dotted line). The measurement 
error has little effect on the threshold value. 

Ft at different probabilities of measurement errors for 
a quantum memory implementing the three qubit bit- 
fiip code. Clearly, Ft is insensitive to the measurement 
errors even when the probability of measurement errors 
is significantly higher than the noise strength 71. The 
probability of the measurement error as high as 5% has 
only minor effect on the threshold value. 

B. Effect of a collective bath 

A distinct feature of our noise model is the ability to 
describe the effect of a collective bath, in which all qubits 
are coupled to the same environment. Such an envi- 
ronment is relevant in physical implementations such as 
trapped-ion quantum computers, where qubits are cou- 
pled to the same collective phonon modes 0, . The 
effect of a collective bath on the fault-tolerant QEC is 
an interesting topic. Because a collective bath seems 
to contradict the idea of uncorrelated and stochastic er- 
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rors that is the foundation of fault-tolerant QEC, sev- 
eral authors have suggested that collective decoherence 
has to be avoided for fault-tolerant quantum computing 
Also, in a collective bath the effects of noise on 
different qubits add coherently; as a result, superdeco- 
herence states exist, and might affect the efficiency of 
fault-tolerant QEC 0\. 

To answer this question, we simulate the fault-tolerant 
QEC circuit for the three qubit bit-flip code using a 
noise model in which all qubits are coupled to a com- 
mon bath. The following forms of correlation functions 
for the stochastic process are used: 

(fc.(t)) - (6Mt)) -0, 

{6e.,{t)6e,{t')) - jo-S{t-t'), 

{5J,{t)5J,{t')) = ^i-6{t-t'), ^ ' 

{6e,{t)SJj{t')) = 0, 

Notice that in Eq. (O, fluctuations on different qubits 
are fully correlated; this reflects the result of coupling to 
a common bath. Figure IHl shows the crash rate constant 
r„ for quantum memories using the three qubit bit-flip 
code with two different types of baths. The crash rate 
curve for the collective bath case is only slightly higher 
than the curve for the localized bath, and there is no 
significant difference between these two lines. This result 
suggests that a collective Markovian bath, which exhibits 
spatial but not temporal correlation, has little effect on 
the efficiency of fault-tolerant QEC. Although superde- 
coherence states do exist when the system is coupled to a 
collective bath, they have little effect on the dynamics of 
the system, because those states represent only a small 
fraction in the whole Hilbert space. The fault-tolerant 
QEC circuit using the five-qubit code was also studied, 
and similar results were obtained. 



C. Repetition protocol 

Our simulation propagates the density matrix of the 
system in the process of computation, therefore, we ob- 
tain the full information about the time evolution of the 
system. By examining the trajectory of the system dur- 
ing the fault-tolerant QEC process, we find the following 
repetition protocol yields the best performance: 

Repetition Protocol B (conditional generation): 

f . Perform the syndrome detection once. If this syn- 
drome is zero, do nothing. 

2. Otherwise, perform the syndrome detection again. 
If the same syndrome is obtained, accept the syn- 
drome and correct data qubits accordingly. 

3. Otherwise, no further action is taken. 

Figure H shows the crash rate constant r„ for quantum 
memories implementing three qubit bit-flip code using 




Noise Strength (y, ) 

FIG. 6; The crash rate per computational step Fn for the 
three qubit bit-flip code as a function of the noise strength. 
Curves for the distinct bath system (solid line) and collective 
bath system (dotted line) are shown. The result for the col- 
lective bath is close to the result for localized baths. This 
result suggests that collective bath has minor efltect on the 
efficiency of fault-tolerant QEC. 

different repetition protocols. Because the majority of 
the measured syndromes will be zero in the weak noise 
regime, protocol B reduces the amount of time required 
for a fault-tolerant QEC step by a factor of two. As a 
result, the crash rate constant per computational step 
F„ decreases by a factor of two when protocol B is used. 
Similar improvements on the fault-tolerant QEC proto- 
col have been suggested by other groups jlO, II, 2^]. 
The idea behind protocol B is that the syndrome detec- 
tion circuit is complicated and generates extra errors on 
the data qubits, therefore minimizing the number of syn- 
drome detection and accepting a syndrome only when 
two consecutive detections agree on the same syndrome 
improve the efficiency of the fault-tolerant QEC proce- 
dure. 



D. Level of parallelism 

An important factor related to the efficiency of a QEC 
circuit is the level of parallelism in the circuit. The syn- 
drome detection circuit shown in Fig. ^is not optimized; 
for example, for a reasonable physical implementation, 
the first two controUed-Z gates might actually be oper- 
ated in parallel to reduce the operation time. Figure |H| 
shows a compressed version of the syndrome detection 
circuit that has increased level of parallelism. Further- 
more, because the interactions used to implement the 
controUed-Z gate commute with each other (Z and ZZ 
commute), the controUed-Z gate can be done in one step: 

controlled-Z = e— ^i^^/4e"(2i+2=)/4 = ^-^^iz^z,-z^-z,)/4 ^ 

This makes it possible to perform all controlled-Z oper- 
ations in a single pulse. Note that this maximal paral- 
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FIG. 7: The crash rate constant per step r„ as a function 
of the noise strength for the two different repetition protocols 
for a quantum memory using three qubit bit-fiip code. Using 
protocol B reduces the crash rate constant by a factor of two. 
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FIG. 9; The crash rate constant per unit time Ft as a function 
of the noise strength for quantum memories using three qubit 
bit-flip code. Curves for three syndrome detection circuits 
dilferent in the level of parallelism are shown. The solid line is 
for the circuit shown in Fig. |3 the dashed line is the increased 
parallelism circuit shown in Fig. |H1 and the dash-dotted line 
is the meiximal parallelism circuit that finishes all controUed- 
Z operations in a single pulse. We see dramatic improvement 
on the noise threshold values when the level of parallelism is 
increased. The result indicates that by increasing the level of 
parallelism, the threshold value can be significantly improved. 



of parallelism is increased, the number of pathways that 
generate uncorrectable errors decreases. We emphasize 
that because our method of estimating noise threshold 
does not assume maximal parallelism, we can access the 
real threshold value that reflects the limitations of each 
individual physical implementations. 



FIG. 8: A circuit implementing the fault-tolerant syndrome 
detection for the three qubit bit-flipping code. 



VI. CONCLUSION 



lelism design does not generally exist for arbitrary phys- 
ical implementations, and is a special case for our choice 
of model interactions (ZZ coupling). 

Figure El shows the crash rate constant per unit time 
for quantum memories implementing three qubit bit- 
flip code. Results for three syndrome detection circuits 
with different level of parallelism are shown. The noise 
thresholds for the original circuit (Fig. increased par- 
allelism circuit ((HJ, and maximal parallelism circuit are 
about 1.5 X 10-2, 2.3 x 10"^ and 4.6 x lO^^, respec- 
tively. The results indicate that by increasing the level 
of parallelism, the noise threshold can be significantly im- 
proved. Note that the reduction of the operation time in 
higher level of parallelism can not account for all of the 
improvement on the threshold values; because the crash 
rate constant per unit time Ft has been scaled by the 
amount of time needed to complete a fault-tolerant QEC 
step (Ft = F„/r), any difference in Fj is from sources 
other than difference in r. The improvement on the 
threshold value is due to the effect that when the level 



We have applied a noise model based on a general- 
ized effective Hamiltonian to study the effect of noise 
on the performance of fault-tolerant QEC circuits. The 
model includes realistic physical interactions for the im- 
plementations of quantum gates, and describes the effect 
of system-bath interactions by including stochastic fluc- 
tuating terms in the system Hamiltonian. As a result, 
this method simulates quantum circuits under physical 
device conditions, and gives us a full description of the 
dissipative dynamics of the quantum computer. 

Fault-tolerant QEC circuits implementing either the 
three qubit bit-flip code or the five-qubit code were in- 
vestigated, and the noise threshold for quantum memory 
and logical X gate were calculated by comparing the log- 
ical crash rate to the error rate of a bare physical qubit. 
The noise threshold of quantum memories using the three 
qubit bit-flip code and five qubit code is about 2 x 10^^ 
and 5 x 10"'', respectively. The noise threshold of logical 
X gates using the three qubit bit-flip code and five qubit 
code is about 1 x IQ-"^ and 4 x 10"^, respectively. Note 
that in our dimensionless system, these noise strength 
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values should be interpreted as the error rate per unit 
time scale At = 1/e, where e is the strength of the con- 
trol fields. These threshold values are obtained from an 
uniform noise model where magnitudes of storage errors 
and gate errors are the same. This result indicates that 
fault-tolerant quantum computing is possible in systems 
with strong storage errors. A possible scenario for such 
system is the linear nearest-neighbor architecture, where 
only nearest-neighbor interactions arc available for two- 
qubit gates, and excess amount of quantum swap gates 
have to be added to the circuit to perform two-qubit gates 
between qubits distant in space. 

We have also carried out a systematic study on several 
variables that can affect the performance of the fault- 
tolerant QEC procedure for the three qubit bit-flip code. 
Our results show that both collective bath and imper- 
fect projective measurement have minor effects on the 
threshold value. However, the repetition protocol and 
level of parallelism can significantly change the perfor- 
mance of the fault-tolerant QEC procedure. Our density 
matrix results indicate that accepting a syndrome only 
when two consecutive syndrome detections agree (proto- 
col B), which reduces the number of required syndrome 
detection steps, is the optimal repetition protocol. Com- 
pared to the simple majority vote algorithm (protocol 



A), protocol B increases the efficiency of fault-tolerant 
QEC at least by a factor of two. Regarding the level of 
parallelism in the syndrome detection circuit, in general, 
a higher level of parallelism results in a more efficient 
fault-tolerant QEC circuit. The improvement can not 
be fully explained by the shorter operational time for a 
more parallelized circuit; we suggest the major contri- 
bution for the improvement comes from the reduction of 
possible pathways for error propagation. Since the level 
of parallelism is actually limited by available physical re- 
sources in reality, it will be interesting to examine and 
simulate this factor according to a specific physical im- 
plementation of quantum computer (such as ion-trap or 
NMR). 

Finally, we emphasize that without specifying the spe- 
cific noise model and physical device conditions, noise 
threshold values are of little usefulness. Our noise model 
is based on well defined parameters that reflect realis- 
tic device conditions, and provides a full description for 
the dissipative dynamics of the quantum computer. As 
a result, this noise model enables us to access the real 
performance of fault-tolerant QEC for individual physi- 
cal implementations. We believe that such information 
can be useful for the design and optimization of quantum 
computers. 
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